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One-dimensional thermodynamic instabilities are phase transitions not prohibited by Landau's 
argument, because the energy of the domain wall (DW) which separates the two phases is infinite. 
Whether they actually occur in a given system of particles must be demonstrated on a case-by- 
case basis by examining the (non-) analyticity properties of the corresponding transfer integral (TI) 
equation. The present note deals with the generic Peyrard-Bishop model of DNA denaturation. In 
the absence of exact statements about the spectrum of the singular TI equation, I use Gauss-Hermite 
quadratures to achieve a single-parameter-controlled approach to rounding effects; this allows me 
to employ finite-size scaling concepts in order to demonstrate that a phase transition occurs and to 
derive the critical exponents. 
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The absence of phase transitions in one-dimensional 
systems is generally understood in terms of Landau's ar- 
gument pj , according to which, macroscopic phase coex- 
istence - and, by implication, a phase transition - can- 
not occur because the system splits into a macroscopic 
number of domain walls (DW); the spontaneous split is 
favored by entropy, which more than compensates for the 
energy needed to create the DWs. 

Landau's argument provides us with a guide to excep- 
tions from the general rule. For example, in lattice sys- 
tems with long-range harmonic interactions of the Kac- 
Baker type and a 4 on-site potential, where a phase 
transition does occur at a finite temperature 0] , the DW 
energy diverges, and therefore Landau's "no go" argu- 
ment is not applicable. A similar situation arises in the 
generic instability model described by the Hamiltonian 



H = 



E 



- 2 vl 



JftiVn ~ Vn-lf + V(y n ) 



(1) 



where y n and p n are the transverse displacement and 
momentum, respectively, of the n-th particle, V(y) — 
(e~ y — l) 2 is an on-site Morse potential and R is a pa- 
rameter which describes the relative strength of on-site 
and elastic interactions; all quantities are dimensionless. 
The model has been proposed in a variety of physical con- 
texts, such as the wetting of interfaces Q and the ther- 
mal denaturation of DNA0|. In the case of the Hamil- 
tonian , the DW is a static solution of infinite energy 
which interpolates between the stable minimum and the 
metastable flat top of the Morse potential ||. Therefore, 
Landau's argument cannot be invoked to exclude a phase 
transition. Whether a phase transition occurs or not can 
only be definitively decided by an exact calculation of the 
thermodynamic free energy. 

In general, thermodynamic properties of Hamilto- 
nian systems belonging [6j to the class Q can be cal- 
culated exactly by the transfer integral (TI) method. 
Standard texts in statistical mechanics impose restric- 
tions in the type of admissible on-site potentials, e.g. 



\im_ y ^± 00 V(y) oc |y| CT , tr > 0; such a restriction - 
which explicitly excludes - is useful in the sense that 
it represents a sufficient condition for the existence of the 
partition function; at the same time, it enforces the an- 
alyticity of the free energy as a function of temperature, 
and therefore, the absence of phase transitionj^, l2l LLiJ- 
In fact, the crucial step in formulating the TI thermody- 
namics of demands the weaker condition of existence 
of a complete, orthonormal set of eigenstates of the - 
possibly singular - integral equation 

dy> e~^T^'-^ 2 G(y,y')Mv') = &Mv) • (2) 
where, in general, 

G(y,y')=e-i v ^ +v ^y^ , (3) 

and T is the temperature. The limiting case V = 
(harmonic chain), with its continuum, doubly degener- 
ate spectrum of plane waves illustrates the above argu- 
ment. In the more general case of the Morse-like po- 
tentials V(y), Eq. J2J can be shown to be singular be- 
cause the corresponding kernel is, similarly, non-Hilbert- 
Schmidt[T3. I am not aware of a general proof that a 
complete orthonormal set of eigenstates exists for this 
class of Hamiltonians; assuming however for a moment 
that this is the case, a phase transition (instability) sce- 
nario is possible if the spectrum contains a discrete and a 
continuum part and the gap between them continuously 
approaches zero at a certain finite temperature, i.e. the 
longitudinal correlation length £ diverges 01 • This is ex- 
actly what happens if we use the gradient-expansion ap- 
proximation (valid for R << 1 in the temperature range 
1 << T « 1/R [l3^ to map © to a Schrodinger-like 
equation. The validity of such a mapping is certainly 
questionable at large values of R. Therefore, it is legit- 
imate to enquire about independent - and more general 
- methods of deciding whether a phase transition occurs. 
In the absence of exact statements about the spectrum 
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of J5J , previous studies have taken a pragmatic approach 
in the verification of the scenario described above; for ex- 
ample, in Q the integral on the left hand side of was 
cut off at a large positive value of y — y ma x and evalu- 
ated on a grid of a given size. This procedure effectively 
approximates by a real, symmetric, matrix eigenvalue 
problem. The numerical procedure is considered satisfac- 
tory if the results do not depend on two large parame- 
ters: the cutoff and the grid size. Other authors [Tl| have 
applied a Gauss-Legendre quadratures procedure to ap- 
proximate the integral in @; although this is somewhat 
more efficient from the numerical point of view, it still 
leaves two large parameters to be dealt with. Therefore, 
the nature of the approach of the matrix eigenvalue prob- 
lem to the limiting singular equation (J2J) remains some- 
what obscure; as a result, the skeptic may ask|l5|: does 
a phase transition really occur in the system defined by 
the Hamiltonian (Q? 

In the present note, I exploit the presence of the Gaus- 
sian factors in the kernel, and approximate the integral 
in the left-hand-side of © by using a Gauss-Hermite grid 
of size N, i.e. 



dye v f(y) 



N 

£ 

n=l 



w„f(y n 



(4) 
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FIG. 2: The magnitude of the gap minimum (circles, right y- 
axis scale) approaches zero as the system size goes to infinity. 
The sequence of the temperatures corresponding to the gap 
minima, T m (L) (diamonds, left y-axis scale), can be used to 
provide an estimate of the critical point T c . 
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FIG. 1: The gap between the two lowest eigenvalues of the 
matrix eigenvalue problem 0, for a variety of N values. For 
a given N, the gap has a minimum at a certain temperature 

where positions and weights are given by the appropriate 
Gauss-Hermite quadratures routine. Besides the obvious 
advantage of eliminating the cutoff from the numerical 
integration, this allows me to identify the largest of the 
Gauss-Hermite roots, jjn ~ {2N + l) 1 / 2 = L with the 



"transverse size of the system" and employ finite-size 
scaling concepts. In this fashion, the singular integral 
equation is approximated as the N — > oo limit of the 
sequence of N x N matrix equations. 

I use "rescaled" variables, i.e. y = py, p — (2i?T) 1 / 2 , 
and divide both sides of (J2J by py/n. This transforms (0) 
to the matrix form 



N 

EM 



e,/T A , 



(5) 



where 



Di 



1/2 e y'v 3 e -(y'-y 3 r 



G{ P y\pf) (6) 



and A u /(2ttRT) 1 / 2 = e-^' T . The advantage of the lat- 
ter rescaling is that the "harmonic background" of the 
free energy has now been absorbed in the prefactor; the 
lowest of the e„ 's expresses the nontrivial part of the free 
energy. 

I have solved numerically |lq | the matrix eigenvalue 
problem © for R = 10.1 @, N = 256, 384, 512, 2048 and 
temperatures in the range 0.85 < T < 1.30. Results 
for the difference between the two lowest eigenvalues are 
shown in Fig. ^ At any given size L, the gap has a 
minimum Ae m (L) at a certain temperature T m (L). Fig. 
12 demonstrates that (i) the value of the gap approaches 
zero quadratically as L — > oo to within 10~ 5 and (ii) 
the sequence of T m (L)'s also approaches a limiting value 
T c = 1.2276 quadratically. 
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FIG. 3: Finite-size scaling of the order parameter. The scaling 
variable reflects the choice v± = 1. Data from the neighbor- 
hood of T c (i.e. \t\ < 0.1 for different values of N tend to fall 
on the same curve; this confirms the choice of vj_ ; the asymp- 
totic slope (dotted line) reflects the property vx = —ft = 1 
(cf. text). 



FIG. 4: "Reduced" fluctuations of the order parameter for 
three different values of N. The common intersection pro- 
vides a further method to estimate T c from finite-iV runs. 
Details of the intersection are shown in the inset, along with 
the estimate T c = 1.2276 (dashed line) obtained above (cf. 
Fig.©. 



I identify the limiting temperature T c , where the spec- 
tral gap of the limiting, infinite-dimensional matrix eigen- 
value equation JSJ vanishes, as the transition tempera- 
ture of the original TI equation (J2J. 

Near the critical temperature T c , the various thermal 
properties of a finite-size system exhibit the competition 
of two transverse length scales: the size of the system 
L (here: L = (2N + l) 1 / 2 ) and the transverse correla- 



1 1/2 



(y 2 ) (yf 



1/2 



If © 



tion length £l = [((Sy) 2 )] 

has the same critical properties as the Schrodinger-like 
equation derived from it within the gradient-expansion 
approximation (e.g. we expect, in the limit of infi- 
nite L, a transverse correlation length £j_ cx \t\ v± and 
an order parameter (OP) (y) cx |t|^ with /3 — — v± = — 1 



and t = T/T c - I. 
system scales as 



Then the order parameter in the finite 



:.'/>l = £/i(|- 



(7) 



where /i(0) = const , and fi(x) cx l/x if x >> 1; the 
first property follows from the requirement of bounded, 
nonzero OP at t = CP and finite L, and the second 
from the requirement of an L-independent limit at values 
L >> the second property guarantees that (3 = — 
as expected. Fig. [3] shows that numerical results ob- 
tained for three different values of N scale properly if u± 
is chosen to be equal to unity. 



Similarly, the OP fluctuations scale according to 



L 



(8) 



where now /a(0) = const, f2(x) cx l/x if x » 1 (cf. 
above, following Eq,EJ|. 

As a consequence of Eqs. J7} and ©, the ratio 



(Sy) L 



— const. 



(9) 



at t = and any L. This provides a convenient graphical 
rule for locating the critical point (cf. Fig. 0] and Ref. 
H3)- The rule is valid as long as (3 = —v_\_, i.e. for 
both 2nd and 1st order instabilities - in the latter case of 
course only those with a continuously divergent OP . 

The finite-size scaling of the gap is described according 
to the Ansatz 



(10) 



Ae L (t) = L-'f G ( — 



where now /g(Q) = const, /g(x) cx x 2 if x » 1 (cf. 
above), and as a result, Aeoo(i) oc with v = 2. 

Numerical results shown in Fig. [S] demonstrate the 
validity of the Ansatz (|10|) . 

I conclude by presenting some typical results of the 
spectra of (J5J at temperatures above T c . Fig. shows 
the values obtained for e n , n = 1,2,3, for T — 1.5 and 
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FIG. 5: Finite-size scaling of the gap. 



various N. For comparison, I have also plotted the cor- 
responding results obtained in the absence of the Morse 
potential (harmonic crystal); note that the system size 
in this case is twice as large, since there is no repulsive 
barrier at negative y. The figure demonstrates that in 
the limit of large L the spectra of both systems behave 
as e n ex (n/L) 2 . In other words, a detailed analysis of 
the spectra of JHJ can be used to demonstrate that the 
thermodynamic properties of the high temperature phase 
coincide exactly with those of the harmonic chain. This 
completes the thermodynamic description of the insta- 
bility of the paricle lattice system as the transition from 
a confined to a deconfined state. 

In summary, I have demonstrated that it is possible to 
view the singular TI thermodynamics of one-dimensional 
lattice systems with a nearest-neighbor harmonic cou- 
pling and a Morse on-site potential as the limit of a se- 
quence of finite matrix eigenvalue problems. The finite- 



size scaling properties of the sequence are consistent with 
the universality hypothesis; in other words, the critical 
exponents of the limiting system with R = 10.1 are all 
identical with those obtained via the gradient expansion 
and the resulting Schrodinger-like equation (under the 
condition R « 1). The procedure described - and, in 
particular the vanishing of the gap in the limit of infinite 
system size - constitutes in effect a "proof that a phase 
transition occurs within the framework of the exact TI 
thermodynamics. 

I thank M. Peyrard, T. Dauxois and J. Jackie for help- 
ful discussions and comments. 
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FIG. 6: The three lowest eigenvalues of © at T = 1.5 (full 
symbols) for various N. For comparison, the three lowest 
eigenvalues of JjjJ with V = (harmonic chain) are shown 
(open symbols). The onset shows the same data plotted as 
e n /{n/L) 2 vs. 1/L; it demonstrates numerically that both 
sets of eigenvalues behave as (n/L) 2 at large L. 
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